General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



A D-AOl 1 098 

CO OBSERVATIONS OF THE EXPANDING ENVELOPE 
OF IRC+10216 

T. b. H. Kuiper, ct al 

Owens Valley Radio jservatory 


Prepared for: 

Office of Naval Research 
National Science Foundation 

National Aeronautics and Space Administration 


1975 


DISTRIBUTED BY: 


National Technical Information Service 
U. S. DEPARTMENT OF COMMERCE 




ADA011098 


\yBQo2 


A Preprint 
from the 


. OWENS VALLEY RADIO OBSERVATORY 

« 

(.'alifornia Institute of Technology 
Pasadena, (California 

1 975 


S. CO OBSnRVATIONS OF thH EXPANDING ENVELOPE OF lRC+10216 


by 


T.B.IL Kiiiper 
Jet Propulsion Laboratory 
California Institute of Teehnology 
Pasadena. California 


Ci.R. Knapp and .S.L. Knapp 
Owens Valley Radio Obser\'atory 

and 

Robert L. Brown 

National Radio Astronomy Observatory* 
Cireen Bank. West V^iriiinia 


NATIONAL TECHNICAL 
INFORMATION SERVICE 

U S -I 

•“Operated by Associated Universities. Inc., under contract 
with the National .Science Fouiulation 

fP 


* *•*» •«« »>«»» l» » 

UuCUMENr CONT.'^OL DATA • K ^ 0 

t»< $1 t». nt tmtm»$n j$ n»wf fc«* »r ^ >*•»• •$r^t tit f •/•.<»# !• flrttmt 

T^«~>*»'*^O^rrT»"cu>t*"7T~ CL* .^IHC • HO-» 


■<*«'%••• ii^c *ci««>*vii •w'As#/ 

jOwens Valley Radio Observatory 
California Institute of Technology 
Pasadena, California 91125 


It 


-UN CLASS 


•• • f . f t ^ 

CO OBSERVATIONS OF THE EXPANDING ENVELOPE OF IRC-H0216 






» Au -MOi* «i ff. #•» aii'Nff* 

T.B.H. Kuiper, G. R. Knapp, S. L. Knapp and Robert L. Brown 


• -4t k'OM t 3 • 7 V 

70. TCf^W NO. NACKt 

»6 NO o ^ 

19 

•• c «:•«*«• *c r 3«« T MO 

NOOO 14-6 7-A-009 4-0019 

^ •••OJtlC T MO 
r. 

4 

#8, 1975 

j 

• ft. OTHfO MCNOOr NOttI (Af*r rfhmr ftumff l^ 0 t hm 

• i>i« r MtMu **Om %r*rvMrMr 

This document ha" been approved for public release and sale; its distributi 
is unlimited. 

1 r. M« M ? *•« v Movrt 

ta. tooNtootNO MiLiranvACTivirT 

Office of Naval Research 


We have observed high-sensitivity emission profiles from the J = 0«-l 
transitions of ^^C^^O and ^^C^^O towards lRC+10216. It appears that the 
spherically symmetric uniform mass-outflow model propt sed by Morris (1975) 
is necessary to describe the line profiles. The outflow appears to be 
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ABSTRACT 


We have observed high-sonsiti vity emission profiles from 
the J = 0^1 transitions of and towards IRC+10216. 

It appears that the spherically symmetric uniform mass-outflow 
model proposed by Morris (1975) is necessary to describe the 
line profiles. The outflow appears to be slightly accelerated, 
having a velocity of 15 km/sec at the edges of the CO cloud, 
compared with 12 km/sec for the more centrally confined molecules 
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I I?JTRODUCTIDN’ 

A 

The infrarod source IRC+1021C (Ncuuebauer and Leighton 1969) 
consists of a small (0'.'4) optically thic)< component at 600 K 
surrounded by a larger (2") optically thin component whose temperature 
is 375 K (Becklin o t al . 1969; Toombs ct al . 1972). Spectra of 
the source suggest that it is a late-type carbon star, surrounded 
by obscuring dust presumably produced in the atmosphere of the 
star itself (Miller 1970; Herbig and Zappala 1970). 

The infrared source is surrounded by an envelope producing 
millimeter-wavelength from many interstellar molecules. Emission 
linos from CO (Solomon et al . 1971), CN and CS (R.W. Wilson et a 1 . 
1971; Turner et al . 1973) HCN (Morris et al . 1971; W.J. Wilson 
ct al . 1973), C^H (Tucker et al . 1974), and SiS, SiO and HC^N 
(Morris et al . 1975) have all been observed. The large velocity 
width (24 kn/sic) seen in these linos strongly suggest?; expansion 
of the envelope. 

The size of the cloud omitting these molecular lines has be«n 
measured as 4 40" for HCN (Wilson et al . 1973) and -'211 for CO 
(Wilson et al . 1973; Ulich 1975). Morris (1975) has propc>sed 
a model for the I.RC+10216 molecular envelope in which the molecules 
are contained in an expanding envelope surrounding the star and 
are excited by absorption of infrared radiation into the excited 
vibrational states, the subsequent decay leaving the molecules 
in higher-lying rotational states in the ground vibrational state. 
Those considerations decrease the required envelope mass from 
>» 1 collisional excitation dominates (Morris et al . 1975) 

to about 10 ^^^‘^^(Morris 1975) . 
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As a result of the TR excitation, the ap^^ar^nt siz*. of the 

molecule-omitting regior depends on the transition observed; the 

size of the region increases as J (the quantum number of the upper 

rotational state) decreases and as the permanent dipole 

o 

moment of the molecule, decreases (see the discussion by Morris, 

1975) . Thus the CO-emitting region is very much larger than that 

of any of the other observed molecules. In particular, it is 

resolvable by the 1-arcmin beam of the NRAO 11-meter telescope, 

so that the CO line shapes should show the effects of partial 

resolution discussed by Morris (x975). 

In the present paper, we report high-sensitivity observations 
12 13 

of CO and CO in IRC+10216, including partial mapping in the 

13 1213 

CO line. The line shapes of CO and CO towards the center 

of the source appear to require a uniformly-expanding envelope 

mode 1 . 

II. THE OBSERVATIONS 

The observations described herein were made in Marc' 1975 
with the NRAO 11-meter telescope at Kitt Peak, Arizona. "^he J = 0*-l 
transitions of ^^CO and ^^CO at 115.2712 ai'd 110.2014 GHz res- 
pectively were measured; the observations reported here were all 
made in clear weather. The telescope beam efficiency was ^^60% 
and the half-power beamwidth was ^^65%. The receiver used the NRAO 
cooled mixer, giving effective single-sideband system temperatures 

(which incorporate factors due to the atmosphere and the beam 

12 13 

efficiency) of ~1200K and ~900K at the CO and CO frequencies 
respectively. The system temperature and atmospheric attenuation 
were calibrated frequently using an oscillating vane (see Penzias and 
Burrus 1973) . The observations were made by position switching 


a reference point 
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1 ° from the source in azimuth, and typical integration tines of 

about 1 hour per point were used. The observed position was 

“l950 * 09^45"*16f8, ^1950 “ +13°30’40" for both and ^^CC; 

in addition, points 1 arcmln north and south of this position 

were observed in ^^CO. The linos wore observed with two 256 

channel filter banks, giving resolutions of 0.1 MHz (—0.26 km/soc) 

and 0.25 MHz (■^O.SS km/soc) per channel respectively. 

The observations are summarized in Figures 1 through 3 and 

12 13 

Table 1. In Figures 1 and 2 the CO and CO profiles centered 

on IRC+10216 are presented; in Table 1 we list for each line the 

* 

values of T , v and Av. These quantities are defined as follows: 

A C 

* 

T^ is the Rayleigh-Jeans equivalent brightness temperature in 

Kelvins (on this scale the peak temperature of Orion A is 60K - 

see Ulich and Haas 1975) ; v is the velocity with rcsoect to 

c 

the LSR; and Av is the velocity half-power width. In I'igure 3, 
the ^ ^CO profiles measured north and south of the central point 
are presented. 

Ill ANALYSIS OF THE LINE PROFILES 

T!ie observational profiles presented in Figutes 1 and 2 show 
qualitatively good agreement with the model profiles predic.ed 
by Morris (1975) and Mufson and liir.zt (1975) for a spherically 
symmcftric, uniformly expandinc envelope. We apply ”, orris' 
description to our observations because it is more a.nenable to 
analysis. The treatment of Mufson and liiszt offers a more 
intuitive understanding and has been included where appropriate. 

In the following discussion, wo use the basic not.atior. 
introduced by Morris (1975): r = radial distance from the center 



of thn source; n = projection of r on the plane of tho sky; 

7. - projection of r along the line of sight; V = expansion velocity 
of source; v = projection of expansion velocity along the line 
of sight; 6z(p,v^ ) = distance element along the line of sight 

for which the 1 ine-o f-sight velocity falls v/ithin ^»v of v^ ; 

G_(p) = antenna pattern with a halt-power boamwidth of R; T (p,v, ) 

= brightness temperature in the spectral line. 

The following relations will be used, explicitly or implicitly, 
in the subsequent derivations: 


= (z' 


2r- 
? ) 


(la) 


V|^ = Vz/r 
p/r = [1 - (V 

- Sv/f^= ^v£(i 


(lb) 

(l c) 

(l d) 


The optical depth along tho line of sight is defined as: 


T(p,v,) 



N ( r ) dz 


(2a) 


where >(r) is the absorption coefficient, is tho local width 
of the line, and 'J(r) is tho CO density. Since* the only 
contribution to ":{p,v^^ ) ct>mes from a path length olonv-*nt <^'7.(Eq 

(equation Id) , 


r (PfV , ) 


«(r)^?(r) g II - (;y-) 

V V 


(2b) 
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We define an equivalent excitation temperature by 


T (r) = ^ 

X k 


. r 

I.' 


, h’C) ^ 

e TF) 

exc , 




\- 1 

-1 

1 J 



v/here T (r) is the excitation temperature of the transition and 
exc 


bb 


is the temperature of the universal backqroend radiation. Tf>c 


brightness temperature is then expressed as 


( 3 ) 


Tg (p, V, ) = (r) exp - 


N 


J 


«x(r) 

TV- 


N(r) dz 


dz 


J 


r 

= Tj^(r) ^ “ <^xp (- ) ) 




since the only contributions to the integral come from the path length 
element ^z. In Fquation (4), r is determined by p and as given in 
Equation Ic. 

The antenna beam pattern is assumed to be of the form 

G(p) = exp , -41n2 ^ ^ (5 

n 3 V 3 J 

and the observed antenna temperature is 


(V ) = 2r 

assuming that the b 
A . Optic a 1 l y Thick 


r 

jG(p) dp 

c 

eam is centered on 
Case ( C O) 


the source. 


(0) 


In the approach of Mufscn and Liszt (1975), the optically thick 
case may be treated by cotuiidoring ccjual velocity surfaces, whicli 
will be the interiors of cones opt*ninq towards t)ic observer \/itli 



If the conns arc not 
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their apices at the center of the envelope, 
resolved by the telescope, 


)/T(0) = rv(r r.ine)^/ftr^ 

» 1 - /V)2 , 

(see Figure 4). 

As the source becomes resolved, lower velocity cones are truncated, 
thus reducing the line center relative to the edges. In the limit of 
complete resolution, the profiles are flat-topped. 

Following Morris' treacment, the brightness temperature (Equation 
4) is given by 

where r is given by Equation Ic. Then 



(V) 


We see that (v,^ ) is characterised completely by the boamwidth H 
and the characteristic source radius Initially, is assumt-d to In- 

uniform out to radius R^, with the integrand in Equation (7) being 
zero for r-^’R^. Tnen, by a suital)lc transformation of Viwiables, 

r 

I 

T, (v ) = ,1 - exp 

M ;* X 

In Figure 1, we sho:/ a function of this form fitted to the observed 
J = 0-1 profile. Since R^, 


-4 ln2 


4 

li 


- (v^/V) 


Uj 


an 


the source radius, and U, the 
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half-powor beamwidth, arr» each approximately 1 arcrr.in (based 
on uhe brightness distribution data of Ulich 1975) they are taken 
as equal. 

To see to what extent the profile shape depends ca the 
assumed form of lot us consider a qaussian distribution; 

* ’’max > 

where is now tt ' half-intensity radius. Then 




= 2-1^ 11 - (v„/V,^ I 


r 


exp(-ln2 ^ ) oxp(-41n2 — L j-2 j 


r; 


B 


+ (B/2R^)^ (1 - (v./V)^)"^)'- 


(9) 


This line profile does not differ from that for uniform excitation 

temperature by more than 15% over the range 0<v^< 0.95V. By 

comparing ‘"'s given by Fquation (5) with the form given 

in Equation (8), we see that when the excitation temperature 

is tapered instead of abruptly terminated, the edges of tfie 

line profile are enhanced sligJitly, but not sufficiently to 

distinguish with present receiver sensitivity, b’e can conclude that 

the line profile agrees well with the assumptions of optical 

tnickness and uniform expansion (or contraction) and that this 

conclusion is relatively independent of the form of T^(r). Since 

12 

we have shown that the CO line is optically thick, we may derive 

the excitation temoerature of the t'', T , from Equation (3), 

* cxc 

to be 7.5 K . 


rdr 
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B. Optir~liy A'hj.n Case 


In the geometrical approach of Mufson and Liszt, an 
thin, uniformly expanding atmosphere may be treated as a 
concentric shells. For each shell (assuming that it is 
the telescope, (see Figure 4), 


’■a*''., 



dv 


2nr sin0.rd9.t 4nr .t 


optically 
set of 
unresolved 


3 

= 0.5 sine d© = dv^ /2V 


by 


So each unresolved shell has a line profile which is flat from 
-V^y ^ V, and consequently, a uniformly expanding, unresolved 
envelope has flat profiles in optically thin lines. If the 
envelope is partially resolved, the low-v.jlocity contribution 
of the outer shells is not seen by the beam, so that the center of 
the profile will be depressed. In the ultimate resolution limit, 
the profile will consist of spikes at r V. 

To treat this case analytically, we use Rquation (4) in the 
limit when 'c(p,v ) ' 1 . Then 


Tjj(?,v ) -oCT^ ^2 V ■ (V, /V)^)'i 

V.’e consider first the simple assumption that constant, 

-2 

that N(r) = Kr , and that either (or both) of those distributions 
truncates at r = R_. This corresponds to Morris' (1975) model. Thus 

(vj 


T^(v ) = 


'AT (r) 


81n2 


K / p 

I 



(-41n2 



p dp 
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where p^(v ) = 1 1 


(v„ /V)^J 


j. 

7. 


This roducos to 


(■— ^— J- {2 I Ir^ ^ (1 - (vyv)^f-) (10) 

-7 1 / ^ 2 

where r (y) = — /exp(-x^) dx . 

V 

Figure 2 shows a function of this form fitted to the observed 
J = 0<-l, profile for TRC+10216, again taking R^/B * 1. 
Again, the agreement is excellent. 

This profile is relatively sensitive to the r dependence of 
« T N. For example, if we take a gaussian dependence for 7 (r) 

^ “X 

and an inverse square dependence for N(r), we find 


‘ ? / ^2 II - !v„/V)=lrr oxpl-ln2 ^ 1 

B / r 


,exp[-41n2 ) p dp 


a T 

max 


2 /-.i: n2' 
B 


K 

V 


I (B/2R^)^ + 1 


(V /V) 



(11) 


In this case, the peaks are less pronounced; for - B the oeakt 

rise up to only 1.4 times the value at v =0, whereas in the 

previous case, the ratio was tv/o. We can also see intuitively 

that a steeper power law dependence will result in diminished 

peaks at v’ = IV. Examination of Figure 2 sho-.'s that the rodel 

described in Equation (10) fits the observational profile very 

well within the noise limits, so that, to within a reasonable 

approximation, the quantity does follow an inverse square 

12 

dependence. Careful mapping of the source in the CO line should 
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be made to derive T^(r), thus Faparating the temperature and 
absorption coefficient/density distributions. 

C . Profiles Offset from Center 

Figure 3 shows ^ ^CO line profiles observed one arcminute 
north and south of the center of IRC+10216. Although these 
positions do not easily lend themselves to an analytical treatment, 

we can draw somt' -iualitative conclusions from their appearence. 

1 3 

From the central CO profile, as well as on physical 
grounds, we conclude that the line is optically thin eveywhero. 

For a source cylindrical along the line-of-sight with uniform 
density and T , we would expect T.(0) to be about half as strong 
as from the central profile because of the relative sizes 

of the source and the beam. For a spherical source with an 
inverse -square density distribution, should be significantly 

less than half as strong as the center position. In fact, the 
observed values of ^ro slightly more than naif of the 

central value (see Table 1). This suggests that the density 
distribution exceeds an inverse-square law near the outer parts 
of the source. 

We note also the absence of high- tuid low-velocity spikes 
in the off-center profiles, which is in accord wi^;h the model 
predictions. The spikes arise from the lines-of-sight v/hich 
pass near the center of the source, and which receive very little 
v/eiqht when the antenna is positioned almost a full beamwidth 


from the source center. 



12 


D. The Content of tho t= :nvelo{G 

While the present observations cannot verify all the details 
of the model, they are consistent with constant excitation temperature 
and inverse-square law density dependence in the TRC+10216 
molecular envelope. The previously derived expressions lead to 

a K = 0.0847 ~ = 1.3 x lo'^^ cm"^ soc"' (12) 

c 1770177 

where V = 1.5x10^ cm/sec, = 5.5x10^^ cm (assuming that the 

source distance is 290 pc (Herbig and Zappala 1970]), is 

the antenna temperature measured at the peaks, and 

12 

is the C.O antenna temperature at the line center. The total 

13 4 

number of CO molecules in the envelope is then given by 9x10 /' . 

IV DISCUSSION AND CONCLUSIONS 

12 13 

V/e have presented high-sensi tivi ty CO and CO line profiles 
observed towards IRC+10216. The profile shapes are consistent with 

IT 13 

the “CO line being optically thick and the CO line thin; the 

^ ^CO line is double-peaked. These profile shapes are consistent 

with those predicted from a partially-resolved, expanding envelope, 

with constant excitation temperature, constant radial expansion 

velocity, and an inverse-square law density dependence (i.e. constant 

mass loss from the central star) . This is the model developed 

for the molecular envelope of IRC+10216 by .Morris (1975). Our 

12 

analysis shows that the shape of the CO lino is deter.Tiined 
largely by the motion of the gas, and is relatively insensitive 
to the distribution of excitation temperature within thv* envelope. 
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The CO excitation temperature in the envelope is 7.5 X. 

Because the molecules are excited by infrared radiation 
from the central star, the CO envelope is much larger than that of 
any of the other observed molecules (e.g. HCN, SiO, etc.). Our 
observations show that the expansion velocity of the CO is 15 
km/sec, while Morris (1975) and Morris et al . (1975) observed 
a velocity of 12 km/sec for the other molecules. These results 
suggest that (the magnitude of) the radial velocity increases 
slightly with distance from the central star; accelerated 
expansion due to radiation pressure on the dust is a possible 
explanation. The flow in the envelope is highly supersonic. 

The shape of the optically thin ^^CO profile is consistent 
with an inverse square law dependence of the density. The inten- 
sities of the ^^CO lines measured one beamwidth away from the 
source center suggest density enhancement in the outer parts of 
the envelope, perhaps where the outflowing gas encounters the 
interstellar medium. 

In all of the above discussion, we have tacitly assumed 
that the observed radial velocity flow is due to expansion of 
the envelope; but the observations are incapable of distinguishing 
between uniform expansion and contraction. In the case of 
IRC+10216, however, we can infer from other considerations (such 
as the observation of the 4.'^^ CO absorption lines by Geballe 

I 

et al . [1973] ) that the envelope is in fact expanding. 

Finally, in order to determine the CO mass of the envelope, 

12 13 

and the C/ C isotope ratio, the whole cloud should be mapped to 

12 13 

high sensitivity in both the CO and CO lines. Such observations 
should also be valuable for studying the temperature and density 
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run in the envelope, and the intoractien between t lie outt'lowinq 
gas and the interstellar mediiin. 
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TABLE 1: Observed Parameters for the CO Lines in IRC+19216 


* 


Line 

T^(peak) 

(km/sec) 

Av (km/sec) 


4.2 

-26. Oil. 0 

31.01 1.0 


0.2 4‘^t.02 

-26. Oil. 0 

29. Oil .0 

^^CO-1 'N 

0.14 !:.02 

-2L. 4il.O 

27. 3:1.0 

^^CO-l'S 

0.18 t.03 

-25.4tl.O 

29.3:1.0 



Refers to the central part of the profile at v = -26 km/sec. 
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Figure 1 


Figure 2 


Figure 3 
Figure 4 


FIGURE CAPTIONS 


12 16 

: C O J * 0<-l profile at the central position of IRC+10216 

observed with a resolution of 100 kHz. The frequencies 
refer to the source rest frame at -26 km/sec. The dashed 
line represents a uniform mass-outflow model in which the 
excitation temperature is assumed constant throughout the 
envelope (Equation 8) . 

: J = 0*-l profile at the central position of IRC+10216, 

observed with a resolution of 250 kHz. The frequencies 
refer to the source rest frame at -26 km/sec. The dashed 
line represents the uni form mass-outflow model with inverse- 
square-law density dependence given by Equation 10. 

: ^^C^^O line profiles measured (a) 1' north and (b) 1' south 

of the central position of IRC+10216. The resolution is 250 kHz. 

: Schematic representation of eui expanding, spherically 

symmetric source distribution, after Mufson and Liszt (1975). 
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